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Abstract. We evaluate the spin squeezing dynamics of N independent spin- 1/2 
particles with exchange symmetry. Each spin interacts with its own reservoir, and the 
reservoirs are independently and identical. The spin squeezing parameter is analytically 
calculated with different kinds of decoherence. The spin squeezing property vanishes 
with evolution time under the Markovian decoherence. Whereas coupled to the non- 
Markovian decoherence channels, the spin squeezing property collapses and revives 
with time. As spin squeezing can be regarded as a witness of multipartite entanglement, 
thus our scheme shows the collapse and revival of multipartite entanglement under 
non- Markovian decoherence. 
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1. Introduction 

Entanglement is a kind of quantum correlations that has been playing a central role in 
quantum information and communication theory. Quantum correlations have attracted 
much interest as fascinating demonstration of nonclassical phenomena and they have 
also found many promising applications such as achieving interferometric [U EJ EJ E] and 
enhancing the ratio-to-noise ratio in spectroscopy [5j |6] beyond the standard quantum 
noise limit. The spin squeezed state is one kind of quantum correlated states [7] with 
reduced fluctuations in one of the collective spin components, with possible applications 
in atomic interferometers and high-precision atomic clocks. Entanglement is based 
on the superposition principle combined with the Hilbert space structure, while spin 
squeezing is originated from another fundamental principle of quantum mechanics — the 
uncertainty principle. It is found that spin squeezing is closely related to and implies 
quantum entanglement [31 El El El EH]- As a multipartite entanglement measure spin 
squeezing is relative easy to be operated and measured. 

However, realistic quantum systems are not closed and therefore it is of fundamental 
importance to study the quantum correlations when the system loses its coherence due 
to the information flows from the systems to the environment [XT]. The unidirectional 
flow of information in which the decoherence and noise act consistently characterizes 
a Markovian process. The dynamics of entanglement in open systems was broadly 
studied [12]. A peculiar aspect of the entanglement dynamics is the well known 
"entanglement sudden death" phenomenon [131 El EH] an d recently it is investigated 
the "sudden death" of spin squeezing under a Markovian process [HI [17]. However, 
there are some systems such as soft- and condensed-matter systems strongly coupling 
to the environment and the coupling leads to a different regime where information also 
flows back into the system from the surroundings, which characterizes a non-Markovian 
process. Memory effects caused by the information flowing back to the system during 
a non-Markovian process can temporarily interrupt the monotonic increase or decrease 
of distinguishability such as spin squeezing parameter. In this paper we will study the 
spin squeezing dynamics of N independent spin-1/2 particles with exchange symmetry 
coupled to individual non-Markovian environments. 

2. Spin Squeezing Definitions 

We consider an ensemble of N two-level particles with lower (upper) state \\) (|t))- 
Adopting the nomenclature of spin-1/2 particles, we introduce the total angular 
momentum (i.e., Bloch vector) 



N 




(1) 
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where S{ = ~&i = | (tl ~~ 14) j (41 V At this point, it is convenient to introduce the 
following spin squeezing parameter [5j [18] 



c2 _ N(AJ H± ) 

" (J? ' 

Here, the minimization is over all directions denoted by n±, perpendicular to the mean 
spin direction n — (J)/\(J)\. If £ 2 < 1 is satisfied, the spin squeezing occurs and the 
iV-qubit state is entangled. 

3. One-Axis Twisted Spin Squeezed States 

Now we introduce one kind of spin squeezed states — one axis twisted spin squeezed 
states. Consider an ensemble of N spin-1/2 particles with exchange symmetry and 
assume that its dynamical properties can be described by collective operators J a , 
a = x,y, z. The one-axis twisting Hamiltonian [191 120] is an Ising-type Hamiltonian 

k = E \w> k ) i 1 - °i) ® i 1 - ' ( 3 ) 

i¥* 

which involves all pairwise interactions with coupling constant f(j, k). 

The one-axis twisted spin squeezed state [2TJ E2, E3] is obtained by the evolution 
of the above Hamiltonian 

i 



—J{j, k)ta J z a 



+r, (4) 



|^>=exp(-i£t)|+)® JV = JJexp 

i¥* 

where |+) = (|t) + |4)) /v2- If we choose the evolution time to satisfy f(j, k)t = mix 
with m an integer, the state \ipt) is a product state. If f(j,k)t = (2m + 1) 7r/2, |"0t) 
becomes a graph state. And for < f(j,k)t < ir/2, \ip t ) is a one-axis twisted spin 
squeezed state with spin squeezing parameter £. 

The spin squeezing parameter of the one-axis twisted spin squeezed state with all 
coupling coefficients satisfying f(j, k)t = a takes this form 

l-(AT-l) [^TW-A] /4 
5 cos 2 ^ 2 « ' [) 

where 

A = 1 - cos^ -2 (2a) , B = 4 sin a cos N ~ 2 a. (6) 
The mean spin direction for the one-axis twisted spin squeezed state is 

n = (cos (Na) , sin (-Na) , 0) , (7) 
and the orthogonal direction is then 

n_i_ = (— cos (p sin (—Na) , cos </> cos (Na) , sin 0) . (8) 

The minimum spin squeezing parameter with respect to a is obtained £ oc l/iV 1//3 shown 
in Fig. 1. 



Non-Markovian Dynamics of Spin Squeezing 



4 



a. 
, o 








4 



7 



InN 



Figure 1. Ln-ln plot of the spin squeezing parameter for one- axis twisted spin squeezed 
states £ v.s. the number of the qubits N optimized with respect to a. 



4. Evolution of Spin Squeezing in The Presence of Decoherence 

We consider a single qubit coupled to an environment which is described by a thermal 
reservoir. The evolution of this qubit is governed by a general master equation of 
Lindblad form 



d 
dt 



X 



Hr,X 



where the reference system is 



H r 



(9) 



(10) 



with A is the strength of the external field proved that the optimal spin squeezed states 
that maximize the sensitivity of the Ramsey spectroscopy are eigensolutions of the 
Hamiltonian [23]. Whereas, the incoherent processes are described by the superoperator 
C: 

b 



1-8) [&+<T-X + X°~+°~- ~ 2<5--X<5+] 



2 

b . 
-s [a-a+x 

2c -b 



(11) 



XO--6-+ - 2o- + xcr- 



[2x - 2^X0-*] , 
o 

with a± = {p x ± i&y) /2. For an arbitrary s, b = and c = 7, the generator Eq. (fTTj) 
describes the coupling between the qubit and a decoherence channel equivalent to a 
dephasing channel. For s = 1/2 and b = c = 7, the qubit is coupled to a depolarizing 
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Figure 2. The time evolution of spin squeezing under non-Markovian and Markovian 
dcphasing. (a) The spin squeezing parameter £ for one-axis twisted spin squeezed 
states coupled to non-Markovian dcphaisng channels v.s. the time t optimized with 
respect to coefficient a and fixed N = 10, 7 = 0.01, and 7/0 = 1. (b) Under the 
Markovian dcphasing, the spin squeezing parameter £ v.s. the time t optimized with 
respect to coefficient a and n{t) = e -0 ' 005 *. 




Figure 3. The time evolution of spin squeezing under non-Markovian and Markovian 
depolarizing, (a) The spin squeezing parameter £ for one-axis twisted spin squeezed 
states coupled to non-Markovian depolarizing channels v.s. the time t optimized with 
respect to coefficient a and fixed N = 10, 7 = 0.01 and 770 = 1. (b) Under the 
Markovian depolarizing, the spin squeezing parameter £ v.s. the time t optimized with 
respect to coefficient a and n{t) = e -0 ' 005 *. 



channel. Whereas, for s = 1 and b = 2c = 7, that is coupled to a decay channel (pure 
damping) . 

Equivalently, one can use the resulting completely positive map £ with x — as 
follows: 

3 

£x = S ^ J Vj^jXO-j- (12) 
3=0 

with x a density matrix for a single-qubit state and Yl^oPj = 1- These noise channels 
are of particles interest in quantum information theory. This class contains for example: 
(i) for po = (1 + 3k 2 ) /4 and p 1 = p 2 = p% = (1 — k 2 ) /4 with k = e -7 * the above 
depolarizing channel; (ii) for p = (1 + k 2 ) /2, pi = p 2 = and = (1 — k 2 ) /2 the 
above dephasing channel. Finally, the decay channel is obtained as 

£x = E QX El + E lX Et, (13) 
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Figure 4. The time evolution of spin squeezing under non-Markovian and Markovian 
damping, (a) The spin squeezing parameter £ for one-axis twisted spin squeezed states 
coupled to non-Markovian damping channels v.s. the time t optimized with respect to 
coefficient a and fixed N = 10, 7 = 0.01, i]o = 100. (b) Under the Markovian damping, 
the spin squeezing parameter £ v.s. the time t optimized with respect to coefficient a 
and K(t) = e -°- 005t . 



with the Kraus operators Eq = and E\ 




VT 



K 



.2 





The quantum master equations with the time-local structures are also very useful 
for the description of non-Markovian processes. Suppose we have a time-local master 
equation of the form 
d 



dt X 



H r ,X 



+ lC(t) X , (14) 



where JC(t) is a time-dependent generator and takes the similar form in Eq. fTTTT) with 
time- dependent parameter j(t). 

If the relaxation rate j(t) is positive function the generator K{t) is in the Lindblad 
form for each fixed t > 0. In the Markovian regime, information of the qubit state leaks 
to its surroundings and \n{t)\ is a monotonically decreasing function of times. In the 
non-Markovian environment, in contrast, information also flows back into the system 
of qubit state and a revival of distinguishability (here the spin squeezing parameter) 
can be observed in the time evolution. With the definition of non-Markovianity [26J, 
we see that an increase of \n{t)\ leads to a negative rate j(t) in the generator K.(t). As 
an example, we consider the case of a Lorentzian reservoir spectral density which is on 
the resonance with the spin qubit transition frequency and leads to an exponential two 
point correlation function 

/(*) = ^o7e" 7t - (15) 
The rate n(t) is defined as the solution of the integrodifferential equation 

^-K(t) = - / dt'f(t - 0«(0 (16) 



dt 

corresponding to an initial condition k(0) = 1. The parameter 7, defining the spectral 
width of the coupling between the qubit and reservoir, is connected to the reservoir 
correlation time r ~ 7 _1 . For a weak coupling regime r] < 7/2, the relaxation 
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time is greater than the reservoir correlation time and the behavior of n(t) is a 
Markovian exponential decay. In the strong coupling regime T]q > 7/2, the reservoir 
correlation time is greater than the relaxation time and non-Markovian effects become 
relevant [IlEellZTlEHlEaEOlEIlBlEaEl. Thus we obtain 

(17) 



n{t) = e 



dt\ 7 ( dt 
I + - sin I 



where d = ^/2r/ 7 — 7 2 - 

We would be interested in the effect of decoherence on the spin squeezing properties 
of a system including N two-level particles. A decoherence model individual coupling of 
each of the qubits to a thermal bath is considered in this paper, where the evolution of 
the kth qubit is described by the map £k with Pauli operators frj (j = 0, 1, 2, 3) acting 
on qubit k. We are interested in the dynamical evolution of a given one-axis twisted 
spin squeezed state ip of TV qubit in the presence of decoherence. The initial state ip 
evolves in time to a mixed state p(t) given by 

p(t)=£ 1 £ 2 ...£ N \^ t )(^ t \. (18) 

For one-axis twisted spin squeezed states, we consider the three kinds of decoherent 
channels. The modified mean spin direction and the orthogonal direction are calculated 

as 

n' = (cos (At - Na) , sin (At - Na) , 0) , (19) 
n' ± = (— cos (j) sin (At — Not) , cos </> cos (At — Na) , sin 0) . 

Coupling to the individual dephasing channel which is the main type of decoherence 
for a spin ensemble, the spin squeezing parameter of the one- axis twisted spin squeezed 
state evolves to 

^ hit) = cos 2 ^- 2 «' (20) 



1 , / . A 2 \ 1 £? 2 



where 

c = 1 + (t) [n - 1) 1 a - - = ) - > (t) (n - 1) - = m ) 

Then the spin squeezing parameter of the one-axis twisted spin squeezed state evolves 
to 



k z (t) cos ZJV z a 



under the individual depolarizing. Whereas, the spin squeezing parameter of the one- 
axis twisted spin squeezed state evolves to 



!«(*) cos^- 1 ^ - [1 -«(*)] y 



in the individual damping channel. 

Figs. 2-4 show the time evolution of spin squeezing of 10-particle prepared initially 
in one-axis twisted spin squeezed state coupled to individual dephasing, depolarizing 
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and damping channels. We compare the evolution of spin squeezing under Markovian 
and non-Markovian decoherence channels. Here we consider a Lorentzian reservoir 
and thus the deocherence function n(t) can be written as exponential decay e~ 7i / 2 
modified by a periodical time-dependent function cos(dt/2) + 7 / d sm(dt / '2) . In short 
time regime, the spin squeezing property collapses and revives under either of three 
kinds of non-Markovian decoherence. With time increasing, the part of the exponential 
decay becomes more important and the spin squeezing property vanishes finally as that 
under the Markovian decoherence does. 

5. Conclusion 

In summary, we study the dynamics processes of the spin squeezing in a spin ensemble 
in which each spin is coupled to independent identical decoherence channel. We 
analytically calculate the dynamics of the spin squeezing parameters under three 
different types of decherence. As we know the Heisenberg scaling 1/N in the 
decoherence-free case can be achieved. In the presence of Markovian decoherence the 
spin squeezing property of one-axis twisted states vanishes with the evolution time. 
Whereas, in the presence of non-Markovian decohernce and in the short time limit, 
the spin squeezing property collapses and revives with the evolution time due to short- 
time memory effect during non-Markovian processing. With time increasing, the spin 
squeezing vanishes finally even under the non-Markovian decoherence. As spin squeezing 
can be regarded as a measure/ witness of multiqubit entanglement, thus our scheme for 
the first time shows the collapse and revival of multiqubit entanglement under non- 
Markovian decoherence. 
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